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1.  Introduction 

In  a  recent  paper  the  author  [l]  used  a  variant  of 
Frledrichs'  [2]  technique  of  stretching  In  order  to  study 
Irrotatlonal  and  rotational  solitary  waves  In  a  long  straight 
channel  with  arbitrary  cross  section.   The  object  of  the 
present  paper  is  to  extend  the  ideas  of  [1]  so  that  we  can  con- 
sider Irrotatlonal  cnoidal  waves  in  a  gravitating  liquid  con- 
fined to  a  toroidal  channel:   a  channel  whose  wall  is  formed 
by  the  rotation  of  an  arc  around  a  fixed  axis  as  shown  in 
Fig.  1  in  Section  2.   We  suppose  that  the  liquid  is  subject  to 
a  body  force  directed  toward  this  axis.   One  of  the  reasons  for 
considering  an  analysis  of  cnoidal  waves  in  such  a  fluid  body 
is  that  it  may  give  some  insight  with  respect  to  nonlinear 
waves  in  a  liquid  or  gaseous  ocean  confined  to  a  channel  along 
a  parallel  of  latitude  on  a  sphere.   It  is  not  a  priori  clear 
that  cnoidal  or  solitary  waves  exist  in  a  toroidal  ocean  since 
their  existence  in  straight  channels  Is  strongly  coupled  with 
the  fact  that  waves  of  very  large  wave  length  exist  in  such 
channels.   The  infinite  length  of  the  straight  channel,  in  which 
it  has  been  proved  that  cnoidal  waves  can  exist,  stands  in 
contrast  with  the  finite  length  of  the  liquid  filled  toroidal 
channel  which  is  the  subject  of  the  problem  treated  here. 

Section  2  contains  a  formulation  of  our  problem  and  we 
present  there  the  equations  which  are  produced  when  a  stretching 
technique  is  applied  to  the  hydrodynamlcal  equations  for  the 
motion  of  an  incompressible  liquid  in  a  toroidal  channel  when 


the  liquid  Is  subject  to  an  axlally  symmetric  gravitational 
force  directed  toward  the  axis  of  rotation  of  the  channel. 

In  Section  3  we  examine  some  phases  of  the  linear  theory 
of  motion  In  the  channel  which  theory  Is  based  on  the  assump- 
tion that  the  wave  amplitude  is  small.   This  section  contains 
a  discussion  of  the  speeds  associated  with  the  characteristic 
sinusoidal  waves  which  are  admitted  by  the  linear  theory. 

Section  4  is  an  amplification  of  an  analysis  completed 
some  time  ago  in  order  to  answer  certain  questions  raised  by 
J.  J.  Stoker  about  two-dimensional  nonlinear  waves  in  a  body 
of  water  whose  equilibrium  shape  coincides  with  a  hollow  cylinder, 
We  demonstrate  in  this  section  that  a  formal  expansion  can  be 
chosen  such  that  the  first  term,  the  lowest  order  approximation, 
represents  a  cnoldal  wave  profile.   The  first  order  speed 
associated  with  this  wave  is  greater  than  the  speed  of  any  one 
of  the  characteristic  sinusoidal  waves  of  the  linear  theory. 
The  representation  for  the  free  surface  displacement  r\ ,    and 
other  representations  which  we  give  for  the  velocity  components 
and  pressure  stem,  as  it  turns  out,  from  the  classical  equation 

(1.1)  m^Ti^"(cr)  =  m^r\^{<x)r]^{<j)  +m^Ti!^{<r)    . 

This  is  the  well  known  equation  which  arises  in  all  of  the  cases 
noted  above.   Under  certain  conditions  it  has  a  solution  which 
can  be  expressed  in  terms  of  the  en  function,  a  Jacobi  elliptic 
function.   We  find  that  the  coefficients  of  (l.l)  depend  on  the 


solution  ■!/'p(p?C)  of  a  boundary  value  problem  of  elliptic  type. 
If  the  cross  section  of  the  channel  is  not  a  rectangle  i/^{p,ll,) 
depends  on  both  p  and  C  and  hence  the  motion  has  three- 
dimensional  aspects.   We  find  in  every  case  that  the  first  order 
approximation  tj^  to  the  free  surface  admits  no   change  in  height 
in  a  lateral  direction  across  the  channel.   However,  if  the 
cross  section  of  the  channel  is  not  rectangular,  higher  order 
approximations  show  that  the  height  is  not  constant  in  the 
lateral  axial  direction.   Thus  the  theory  presented  here  leads 
to  a  three-dimensional  wave. 

These  results  and  others  are  summarized  and  discussed  in 
more  detail  in  Section  5- 


2.  FormulatJon 

Let  -r(P)  denote  the  vector  directed  from  a  point  P  to  the 
nearest  point  Q  which  lies  on  a  fixed  x^-axis.   If  5  is  the 
density  of  a  particle  at  P,  and  if  r  is  tlie  distance  I'rom  Q  to 
P,  the  vector  -  —  •  x'fP)  defines  a  force  field  which  can  be 
regarded  as  a  gravitational  field.   We  suppose  that  this  field 
compels  an  Incompressible,  inviscid  liquid,  with  constant  density 
5  to  remain  in  a  toroidal  channel.   The  inner  wall  of  the 
channel  is  the  surface  generated  by  rotating  the  piecewise 
smooth  curve 

(2.1)  r  =  a  +  Y(x^)  =  ah  +  Y(x^) 


around  the  x^-axis.   The  letter  a  denotes  a  free  constant  while 
h  denotes  a  unit  of  length.   A  cross  section  of  the  channel 
lying  in  the  XpX  -plane  is  shown  in  Fig.  1,   The  cross  at  the 


(hC-]^,ah 


*■    X- 


Fig.  1 


origin  in  the  diagram  indicates  that  the  x-,-axis  is  directed 

into  the  paginal  surface.   We  assume  that  the  x, ,  x„,  x^-axes 

i    ^    j5 

define  an  inertial  rectangular  reference  frame.   We  suppose  for 
simplicity  that  a  line  parallel  to  the  x-,-axis  intersects  the 
channel  wall  in  no  more  than  two  points;  that  Xp  =  ah  intersects 
the  wall  in  the  points  (hC,,cth)  and  {\\t,    ,o?ci)    such  that 
b  =  Ph  =  h(C2  -  C-,  )  >■  0,  where  b  is  the  breadtl^i  of  the  equi- 
librium free  surface.   We  also  suppose  that  the  maximum  depth 
of  the  liquid  in  the  channel  is  small  compared  to  the  distance 
of  the  bottom  from  the  x-j-axis  wliich  we  assume  is  larce . 
We  will  use  the  half  plane  M  defined  by 


Xg  =  0  ;    x^  >  0 


as  the  basis  of  a  cylindrical  coordinate  system  in  which  a 
point  has  coordinates  (r,9,x^)  where  0  is  the  angle  which  the 
vector  r  makes  with  M.   The  angle  0   is  to  be  measured  from  M 
in  a  counterclockwise  direction  —  the  direction  defined  by  the 
palm  directed  fingers  of  the  right  hand  when  the  thumb  points 
in  the  direction  of  the  positive  x^-axis.   In  terms  of  these 
coordinates  the  equation  of  the  free  surface  of  the  liquid  in 
the  channel  can  be  written  as 


(2.2)  r  =  a+  f (9,x  ,t) 


where  t  stands  for  time. 

Let  the  cylindrical  velocity  components  of  a  liquid  parti- 

6.0  dr 

cle  be  V,  =  -r  -ttt,  the  clockwise  tangential  speed;  v   =  -g-p,  the 

radial  speed;  and  v  =  -r-p^  the  axial  speed.   If  the  gravitational 

force  and  the  pressure  p  are  the  only  forces  acting  on  the  liquid 

the  elementary  theory  of  hydrodynamics  predicts  that  the  motion 

of  the  liquid  is  determined  by  the  continuity  equation 


Sv,    Vp    Bvp    3v- 
F  dd  ~        dr  "Sx- 


-,  ov,    v^    ov„    ov^ 
(2.3)  -  i-^   +  ^  +^  +  ^=  0  ; 


the  momentum  equations 
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the  kinematic  boundary  conditions 


(2.5) 


5Y(x  ) 

^2  =  ^5  -^3r^ 


"3 


(2.6) 


""2--—  ^  ■  ^3^       ^ 


plus  dynamical  conditions  at  the  free  surface;  and  initial 
conditions  at  t  ^  0.   The  vortlclty  vector  is 


(2.7) 


-  "  ^1-1  +  ^2-2  +  ^5^5 


_,^^3 


Svp  •      / ^  Sv    ^^1 

^ir;  ^  ^1  ■^  I  ?  ^F"  +  ^ir:  ^2 
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Here,  k.,  ,  kp  and  k^  are  mutually  orthogonal  unit  vectors  which 

define  a  local  coordinate  system  at  any  point  P.   The  vector  k^ 

E. 
lias   tlie   same   direction  as   the  x-,-axis;   ]£p  =  -^r    is  a  radial 


vector,  and  ]<,  points  clockwise. 


The  above  equations  can  be  written  in  d Lmensionless  form 
if  we  introduce  a  typical  length  in  the  radial  direction  say  h; 
a  dlmensionless  stretching  factor  /e";  and  the  dimensionless 
variables 


cr=  /e 
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^    h 


/^  (f)^/^  t 


Vn 


U   = 
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V  = 
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w  = 


gh 


e  V. 


'gh 


TT 
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Tgh 


a 
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y  =  u 
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In  terms  of  these  variables  the  basic  equations  are 


(2.8) 


e  Su   V   9v   Bw  _  ^ 


e[ 


5u 


u  (3u   1  Stt-,  ,    Su  ,    5u  ,  uv   ^ 
p  d(r   p  60^'  dp  dQ    p 


(2.9) 


u  Sv    u   ,  -,  ,  Stt-i  ,    5v  ,    hv 

—    ^— -  -  +  1  +  ^r —    +  V  -r —  +  W  -r-^ 

p  ocr   p       dp      dp     ac. 
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"-dT 
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u  Sw    Stt-, 
p  30^  +  ^J 


Sw  ,    civ;   „ 

V  ^r—  +  w  ^r^  =  0 

dp      dC 


(2.10) 


V  =  w 


ay  . 


(2.11) 


^t^-^^l^-w|^; 
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V-,  = 


(2.12) 
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When  the  motion  of  the  liquid  In  the  channel  Is  Irrota- 
tlonal  there  exists  a  velocity  potential  ^ip,(T,C)    such  that 
for  a  liquid  particle  with  coordinates  (p,<r,C)  the  velocity 
components  are 


(2.13) 


i 


u  =  - 


V 


h       ,         w  =  (J)^  ; 


and,  by  virtue  of  the  continuity  equation,  ^   must  satisfy 


(2.14) 


^'t'^^   ^ 


crcr 


P   ^pp   ^CC 


The  momentum  equations  In  terms  of  ^   are 


0 


"1  "1     Q 

Vt    2  VVp   ^3  xr      p    p  pp  ^C^Cp 


r 


=  0 


and  these  yield  the  energy  Integral 
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(2.15)       £[h      +  -^  +   p  +7r]  +  -£.  +  -^  =   const. 


This  Integral  requires  the  dynamical  surface  condition 

^  ^  \)  ^ 

(2.16)      e[({)   +  -2-^  +  a+Ti  +  7r]  +  -£-  +  -S.  =   const. 

2p  "^  ^ 

to  be  satisfied  at  the  free  surface  p  =  a  +  rj  (cT, /^,  t)  ,   In 
addition  to  this,  equations  (2.11)  and  (2.13)  show  that  the 
kinematic  condition 


(2.17)  e[Ti^  +  -^]    =   i   -  (j)^Ti 


T  "  ^^J  =  ^p  "  ^f;''c 


must  also  be  satisfied  at  p  =  a  +  rj .   In  terms  of  9,  equation 
(2.10)  is 


(2.18)  (j)   =  (t)^y^ 


which  asserts  that  the  normal  velocity  of  the  liquid  is  zero  at 
the  rigid  channel  wall  p  =  ct  +  y(r). 

The  equations  (2.8)-(2,l8)  are  analogous  to  the  equations 
which  were  recently  used  by  the  author  [l]  in  order  to  study 
solitary  waves  in  a  horizontal  straight  channel  with  arbitrary 
cross  section.   The  presumption  here  is  that  (2.8)-(2.l8)  will 
be  useful  for  the  Investigation  of  nonlinear  waves  in  a  toroidal 
channel.   However,  before  we  examine  this  matter  let  us  consider 
in  the  next  section  some  linear  aspects  of  the  above  equations. 


3.  Linear  Theory.   Irrotatlonal  Motion 

For  the  purposes  of  this  section  no  stretching  Is  needed, 
that  Is,  we  take  e  =  1  so  that  cr  =   9. 

If  the  motion  of  the  liquid  In  the  toroidal  channel  Is 

such  that  u=u,  v=v,  w=w  where 
o'      o'      o 

bu  Su 

^^  =  0  ,    ^TTT^  =  0,    v=0,    w=0 
dT       '69       ^     o      '     o 

then  the  wall  condition  (2.l8)  Is  satisfied  and  we  will  say 
that  the  channel  contains  a  liquid  stream.   As  we  can  see  from 
equations  (2.12)  a  liquid  stream  is  Irrotatlonal  only  if 

const.     c       ^o9 

u  =  =  =  -  . 

o      p 


gh       ^ 

For  this  kind  of  motion  equation  (2.15)  shows  that  the  pressure 
is  a  function  of  p  only  and  (2.l6)  shows  that  the  equation  of 
the  free  surface  is  p  =  a  +  const.   Since  the  constant  a  Is 
arbitrary  the  free  surface  can  be  defined  by  p  =  a  (t]  =  0);  and 
if  we  take  the  free  surface  pressure  to  be  zero  we  have 

^11 

^o  =  "'-p  +lgh  (;t  -  -?)  • 

Let  us  consider  the  problem  of  finding  the  Irrotatlonal 
steady  motion  which  results  when  we  apply  a  slight  pressure 
F(9,C)  to  the  surface  of  a  liquid  stream  as  defined  above.   The 
potential  function  for  the  disturbed  stream  can  be  expressed  as 


10 


while  the  pressure  and  the  free  surface  are  respectively 
represented  by 

2 

TT  :=  TT^  +  TT*  =  a  -  p  +  |-j^  (^  -  -^)  +  TT* 

^   a    p 
and 

p  =  a  +  Tj*  , 

Without  any  stretching,  that  is  witli  e  =  1  In  (2.  l4 )- (2.  l8) , 
the  starred  perturbation  quantities  must  satisfy  the  following 
nonlinear  system.   The  equations 

_2     p    ^pp    ^CC 


and 


(ig)  C(j)g  ^       ((t)p)         ((j)^) 


(5.2)     2 2 +  TT*  +  — £ —  +  — ^ —  =  const. 

2p     p  /gh 


must  hold  at  any  point  occupied  by  a  liquid  particle  while 


p      ygh     ""^   ^  ^ 

and 

2    ^     ^ 
(3.4)  F(0,C)  =  a-  p  +  |-^  (  1     1  )  +Tr^ 

'^    a    p 


must  be  satisfied  at  the  free  surface  p  -   a+r)*,  and 

(3.5)  <t>p  =  i^y^ 

11 


must  be  satisfied  for  p  =  a+y(C)«   These  equations  can  be 
linearized  by  neglecting  products  of  starred  quantities.   If 
we  do  this  we  find  that  (j)  must  satisfy 


(3.6) 


~2^         p    ^pp    ^CC 


for  0  <  0  <  27r  and  (p,C)  in  the  domain  D  shown  in  Fig.  2.   The 


*^  ^ 


Fig.  2 


linearized  free  surface  boundary  conditions  are 


(5.7) 


(3.8)    F{e,c)  = 


-(a\h  -  c^)  ^  +  — % 

a^gh     a  /gh 


+  const, 


and  the  wall  condition  is 


12 


(5.9)  <l>p  =  <t)Jy^  . 

The   elimination   of  r]    from    (5-7)    and    (3-8)    gives 

(3.10)  a(a^gh-  c2)(t)^(a,e,a  +c^<t)9g(ci,e,a    =  a^c/ih  F0(  0,  C)    • 

When  the  pressure  P(0,C)  and  its  derivative  Fg(9,0  are 
small  we  can  expect  on  physical  grounds  that  the  original  non- 
linear system  will  possess  a  unique  solution  for  any  finite  c, 
or,  what  amounts  to  the  same  thing,  for  any  mean  velocity 


c    II   j   dCdp 

D 


dCdp 
D 

of  the  undisturbed  stream.   It  also  seems  reasonable  to  expect 
that  the  solution  of  the  corresponding  linear  system  will 
provide  an  approximation  to  the  solution  of  the  nonlinear'  system. 
It  turns  out,  however,  that  the  latter  expectation  cannot  be 
completely  fulfilled.   We  will  show  that  for  certain  values  of 
c",  which  cannot  be  precluded  on  physical  grounds,  a  solution  of 
the  linear  system  for  a  rectangular  cross  section  does  not  exist; 
and,  in  addition,  for  values  of  c"  in  the  neighborhood  of  such 
values  the  linear  system  can  only  be  satisfied  by  quantities  so 
large  that  the  linearization  is  inadmissible. 

For  the  special  case  of  a  toroidal  channel  with  a  rectang- 
ular cross  section  of  breadth  b  =  Ph  and  depth  h,  the  disturbance 

I  * 
potential  (J)   must  satisfy 


15 


(5.11)  i|£  A.^^p.^-^.o 

for  0  ^  0   <_  27r,    0   <   a-1   <p<a,    -  ^  <    1^    <  ^  and   the  boundary 
conditions 

(3.12)         a(a\h  -c^)(l)p(a,9,0  +c^(t,gg(a,e,C)    =  a^c/ih  Pg(  8,  C  ) 


(3.13) 


<l>p(a-l,0,C)    =   0 


<l)*(p,     ,±  |)    =   0 


Let  us  suppose  that  the  pressure  is  Independent  of  C?  say  given 
by 

(3.1^)  P(0,C)  =  F*{e) 

so  that  we  have  a  two-dimensional  equation  and  boundary  condi- 
tions: 

(3.15)  p^<t)pp(p,0)+  P^p  +  <t>00  =   0 

0   <   9   <   2Tr   ,        a-1   <   p   <   a 


(5.16)        a(a\h-  c^)<t)p(a,0)+  c^(})gg(a,e)    =  a^c/gh  F*(0) 


(3.17)  ^^(a-1,0)    =   0    . 

r 

An  application  of  finite  Fourier  transforms  shows  that  these 
relations  are  satisfied  by 


14 


(3.18)  ^*ip,9) 


mst.  +^  Y. 


«^c/ih[(^)J+  (^)J]  rF*(x)sinj  (x-0)dx 


00 

^    ,  -   .  -TT 

=  cor~ 


j  =  l         (ct  gh  -  c  )sinh  j^  -  jc   cosh  j^ 
where 

(3.19)  c^  =  ^ 


a-l  ' 


and  this  is  the  only  solution  of  the  problem  that  exists 
provided  no  denominator  in  the  sum  is  zero.   If  for  some  j  =  m, 
c  is  such  that 

3      2  2 

(3.20)       (a  gh  -c  )sinh  m^ -  mc   cosh  m^  =  0 

then  a  solution  exists  only  if 


(3.21)  /   F*(x)sln  m(x-0)dx  =  0 

-TT 


When  (3.21)  and  (3. 20)  prevail,  the  solution  is  not  unique  and 
it  has  the  form 

(3.22)    ^*{p,e)    =   const.  +A[(^)"'+  (^)'^]sin  m(e-B) 


TT 


^hi 


00  I 


a'^c/ih[(^)'U  (^)-^]  /   F"(x)sin  j(x-0)dx 


-TT 


"^TT  i^-^  /  2  ,    2v  .  ,  71      T"?    r  TT 

j  =  l         (a  gh  -  c  jsmh  j^  -  jc   cosh  j  t^ 


where  the  prime  in  the  sum  indicates  that  the  m-th  term  is  to 
be  deleted,  and  A  and  B  are  arbitrary  constants.   If  F*(9)  =   0 


15 


and  (3' 20)  holds  the  solution  of  the  homogeneous  system  is 

(3.23)     <^*ip,9)    =   const.  +  A[(^j-)"'  +  («=.l)'^]sin  m(0-B) 

The  values  of   c   which  satisfy  (3'20)  for  m  an  integer 
m  ^  1  ,  that  is, 


(3  24)         c   =  / ^S]l_ 

^-'•'^  ''  ^m    /l  +  m  GOth  m^ 


are  values  for  which   4>   does  not  in  general  exist;  and  further- 
more, the  linearization  of  our  problem  can  be  questioned  if   c 
is  sufficiently  close  to  a   c   because  then  the  magnitude  of 
the  m  '   term  in  the  sum  (3-l8)  is  as  large  as  we  please  if 
(3-21)  does  not  obtain.   Tne  values   c    also  define  the  velo- 
cities of  the  liquid  stream  for  which  in  the  absence  of  surface 
pressure  there  can  be  steady  motions  other  than  the  streaming 
motion,  namely  those  given  by 


* 
o 


+  const.  +  A[(-e-,-)^  f  («--l)"^]sin  m  (S-B) 


c.e 

+  const.  +  RL  ^   1  /   '     \ 

/gh  "-1        P 


The  mean  velocities  of  the  stream  corresponding  to  a  c   are 

m 

given  by 


i6 


m 


a   ,^t2   a^gh  in^  ^ 


^  a-1 


a-1 


1  +  m  coth  m^ 


3i 

e 


o  sinh  -^  •  [c  4-m4  coth  m^  J 


It  Is  interesting  to  notice  that  as  the  radius  of  the  channel 
is  allowed   to  become  infinite,  that  is,  as   a  — >  oo  ,  or 
I  — )   0  ,  we  have,  for  any  finite   m  , 


11m    c   =  '/gh 
m 

Or—*.  00 


which  is  the  so  called  critical  speed  for  uniform  flow  in  a 
straight  horizontal  channel  with  rectangular  cross  section 
of  depth  h. 

For  the  case  of  a  horizontal  channel  with  rectangular 
cross  section  J.  J.  Stoker  [3]  has  shown  that  a  reason  for 
designating  y/gh     a  critical  speed  can  be  found  in  an  analysis 
of  certain  initial  value  problems.   For  example,  if  a  time 
independent  pressure  is  applied  at  time   t  =  0   to  the  surface 
of  a  stream  running  uniformly  at  any  speed   U  in  the  channel 
then,  within  the  framework  of  the  linear  theory,  equilibrium 
initial  conditions  can  be  used  to  find  the  unsteady  distur- 
bance potential.   It  turns  out  that  no  matter  what  the  un- 
disturbed equilibrium  speed   U  is,  the  disturbance  is  finite 
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for  finite   t.   Also,  If  U  7^  /gh   the  disturbed  flow  generally 
tends  to  a  steady  flow  as   t   tends  to  Infinity.   However,  If 
U  -   \/gIi   then  there  can  be  no  steady  state;  the  disturbance 
potential  tends  to  Infinity  at  all  of  Its  space  coordinates 
as   t  approaches  Infinity;  and  this  Is  one  reason  why  \/gh 
Is  called  a  critical  speed.   Stoker  [3]  notes  that  even  if  a 
dlsslpative  force  is  Introduced  so  as  to  cause  the  unsteady 
motion  to  die  out  the  resulting  steady  motion  still  becomes 
unbounded  if   U  =  \/gh.   One  is  thus  led  to  the  conclusion 
that  the  linear  theory  falls  when  the  speed  U  is  /gh  and 
for  this  speed  a  nonlinear  theory  is  requlredo 

In  a  similar  way  it  can  be  shown  that  if  a  time  indepen- 
dent pressure  is  applied  to  the  surface  of  a  liquid  stream 
running  as  defined  above  in  a  toroidal  channel,  then  when 
c  =  c    the  unsteady  disturbance  potential  tends  to  infinity 
for  all  Q     and  p   as  the  time  tends  to  infinity. 

The  values   c   must  satisfy 
m  '' 

c^  -     °^  gh      ^    g^gh  g,     <  g^gh  C 
m    1  +  m  coth  m^  —  £,  +  E,  coth  ?,    ^  +  1 

and  when   a   is  large   a'^gh  |/(^  +  ^  coth  4)   is  close  to 
ot  gh  ^/{^  +  1).   Note  that  the  speed  given  by 


1    1  +  coth  ^ 


c2  ^    a^P;h 
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In  the  speed  of  the  sinusoidal  waves  of  greatest  wave  length 
namely  2iro,.      It  Is  also  interesting  to  observe  that  the  solu- 
tions (3«l8)  for  <^*      is  unique  if 

(3.25)  c^>^f^        . 

When  (3*25)  is  satisfied  the  linear  theory  predicts  that  in 
the  absence  of  surface  pressure  no  characteristic  sinusoidal 
surface  waves  exist. 

4.   Cnoidal  Waves 

Proofs  of  the  mathematical  existence  of  solitary  and 
cnoidal  waves  appeared  only  recently  after  many  scientists 
had  analyzed  these  waves  in  numerous  papers.   Friedrichs  and 
Hyers   [4]  (1954)  and  Lavrentiev  [5]  (194?)  proved  the  mathe- 
matical existence  of  two-dimensional  irrotational  solitary 
waves.   Littman  [6]  (1957)  proved  the  existence  of  two- 
dimensional  irrotational  cnoidal  waves.   These  proofs  were 
given  about  eighty  years  after  Boussinesq  [7]^  [°] ,    (1871, 
1872)  and  Rayleigh  [9]  (1876)  independently  of  each  other 
presented  the  first  mathematical  analysis  of  the  solitary  wave. 
In  other  words  the  proofs  were  given  about  one  hundred  ten 
years  after  Russell  [lO],  [ll]  (1838,  l844)  wrote  the  first 
scientific  reports  on  the  solitary  wave.   Extensive  references 
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to  the  literature  can  be  found  In  Vol.  IX  of  the  Encyclopedia 
of  Physics  (Springer-Verlag)  which  contains  the  exposition 
Surface  Waves  by  Wehausen  and  Laitone. 

The  existence  proofs  mentioned  above  show  that  if  h 
is  the  depth  of  the  channel  these  cnoidal  and  solitary  waves 
travel  with  speeds  in  the  neighborhood  of   U  =  /gh  which  as 
we  have  already  noted  is  called  a  critical  speed  in  the  linear 
theory.   According  to  this  theory  no  waves  can  exist  on  a 
running  stream  subject  to  zero  surface  pressure  if  the  speed 
of  the  stream  is  equal  to  or  greater  than  /gh  ,  whereas  for 
each  speed  less  than  this  value  a  sinusoidal  wave  is  possible, 
and  its  wavelength  tends  to  infinity  as  the  speed  tends  to 
/gh  from  below.   For  other  characterizations  of  the  critical 
speed  see  Stoker's  book  [3]  and  the  paper  by  Peters  and  Stoker 
[12]. 

Two-dimensional  rotational  solitary  and  cnoidal  waves  have 
been  investigated  by  Hunt  [13]^  Peters  and  Stoker  [l2],  Moiseev 
[l4],  Ter-Krikorov  [15],  Benjamin  [16] ,    and  Shen  [17] .   Ter- 
Krlkorov  [15]  proved  the  existence  of  a  type  of  two-dimensional 
rotational  solitary  wave  by  using  the  method  which  Friedrichs 
and  Hyers  devised  for  the  irrotational  case.   In  a  recent  paper 
Peters  [l]  investigated  irrotational  and  rotational  solitary 
waves  in  a  long  straight  horizontal  channel  with  arbitrary 
cross  section.   Under  all  circumstances  considered  in  the  papers 
cited  above  Peters  [l]  showed  that  the  speeds  of  the  nonlinear 
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waves  are  predicted  by  the  linear  theory  to  be  In  the  neigh- 
borhood of  the  speeds  of  the  long  waves  of  that  theory. 

Coupled  with  what  we  found  In  Section  5^  the  results  Just 
mentioned  suggest  that  Irrotatlonal  cnoldal  waves  may  exist  in 
a  toroidal  channel  with  rectangular  cross  section  and  travel 
with  speeds  In  the  neighborhood  of  the  critical  value 


a\h  i 
^    +   ^  coth  ^ 


which  is  nearly  the  same  as 


(4.1) 


a-^gh  ^ 


when  a   is  large.   In  order  to  examine  this  possibility  and 
the  more  general  one  of  cnoldal  waves  in  a  toroidal  channel  of 
arbitrary  cross  section  we  will  use  the  equations  (2.13)  -  (2.l8) 
which,  as  we  saw  in  Section  2,    are  produced  by  applying  a 
stretching  technique  to  the  basic  hydrodynamlcal  equations. 
The  equations  are 


(4.2) 


u  =  --i 


v  = 


w  = 


(4.3) 


(4.4) 


e<J> 


2     p     pp    CC 


.2p' 


+  p  +  TT 


4)2    (t)2 

+   -^  +   -L 
2     2 
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(4.5) 


2p' 


+  a  +  T] 


<t>2    ,2 


(1^.6)     ^  *^Ti^  =  c^p  -  c^^n^ 


(4.7) 


P  =  ^C'C   • 


Equation  (4.3)  must  hold  at  each  point  {p,o,^)      occupied  by  a 
liquid  particle.   Equations  (4.5)  and  (4.6)  must  be  satisfied 
at  the  free  surface   p  =  a  +  ri(a,0   while  equation  (4.7)  must 
be  satisfied  at  the  wall  of  the  channel  defined  by  p   =  a  +  y{l^) . 
Equation  (4.4)  allows  us  to  compute  the  pressure  once   *   Is 
known. 

Our  object  Is  to  examine  the  consequences  of  the  assumption 
that   *  and  "n   can  be  expanded  In  powers  of   £  ,    that  Is, 


(4.8)      <}'(p,a,0  =  *o  ^  ^*1  "^  ^^*2  "^  ^^S  ^ 


(4.9)     T](a,C)=  ^^2^  +   e^^2  ^   '^^   +  ... 


where   <J>  .   and   t]  .  are  Independent  of   e.   [We  know  that  expan- 
slons  of  this  kind  are  asymptotically  valid  for  two-dimensional 
cnoldal  waves  in  a  horizontal  channel  and  we  can  expect  them  to 
be  valid  in  the  same  way  for  the  case  of  a  toroidal  channel 
provided  a  is  large.]   The  idea  is  to  substitute  these  expan- 
sions in  (4.3),  (^.5),  (4.6),  (4.7) j  and  then  equate  coefficients 
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of  like  powers  of  e.   This  leads  to  a  sequence  of  partial 
differential  systems  for  the  determination  of  (j)  ,^,,...; 
and  ■Hnj^p'**'  *   ^^  need  to  observe  here  that  the  constant 
C  in  (4.5)  must  also  be  replaced  by  an  expansion  in  powers 
of  e.   We  will  write 


(4.10)        C  =  C  +  eC.  +  £^Co+  s^C^  + 
^     ■'  o     1      2      5 


The  zeroth  order  system  comes  from  substituting  e  =  0 
in  (4.3),  (4.5).  (4.6)  and  (4.7).   We  find  that  ^    (p,5',C) 
must  satisfy 

(4.11)         ct)    +-^  +  ^^n  =   0 

opp    P  "^^ 


in  the  domain  D  shown  in  Fig.  2  in  Section  3'   This  domain 
is  bounded  byL:   p  =  ct+y(iC);  and  a  segment  of  the  line 
p  =  a  determined  by  the  points  where  L  intersects  p  =  a. 
We  assume  that  L  is  a  piecewise  smooth  curve  which  crosses 
p  =  a  in  just  two  points  [^    ,a)    and  {Cp,ci-)    such  that 
p  =  C2-  ^1  ^  0  and  C^  -=  0.   The  equations  (4.5)  and  (4.6) 
show  that 

(4.12)         \>^    ia,(J,0  +  ^l^ia,(5,0    =   C^ 


(4.13)  <t)^p(a,(y,C)  =  0 

and  (4.7)  shows  that  the  normal  derivative  i   must  be  zero 
along  L.   Since  the  above  requires  the  normal  derivative  of 
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(j)   to  be  zero  along  both  L  and  S  it  can  be  proved  that 


o 


'i^o  =   ^o^^) 


is  the  only  solution  without  singularities  at  the  boundary. 
Then  from  equation  (4.12)  we  find 


Co  =  0 


In  order  to  fix  F  {<s)    we  must  pass  to  the  next  system. 

The  first  order  system  comes  from  equating  coefficients 
of  e  in  (4.3),  (4.5),  (4.6)  and  (4.7).   The  result  of  doing 
this  shows  that  (()-,  must  satisfy 


(f),  F"(cr) 

(4.14)         ^,    +  _L£  +  ^    =  _  _2_^ 


in  D,  and  for  S  we  obtain  the  conditions 

^^ 
(4.15)  -^  +  a  =  C. 

2a 


(4.16)  (t)^p(a,o-,0  =  0  . 


In  addition  we  find  that  the  normal  derivative  ({>-,   must  be 
zero  along  L.   Equation  (4.15)  shows  that  Fq('^)  is  a  constant 
for  we  have 
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[F;(or)] 


2a 


2 +  a  =  C^ 


(^.17)        [F;(or)]^  =  (C^-  a)2a2  .  ^ 


Consequently  ()),  must  satisfy 


^- 


^Ipp    P    ^ 


=  0 


-PP  "  P  "  ^CC 

in  D  and  since  i-,   is  zero  on  L+S  we  see  that 

^In 

(^.18)         <t^i((r,p,C)  =  F^(cr)  . 

The  second  order  system  shows  that  ({)„  must  satisfy 

(*.i9)         ^p  +  ^  +  4,^^^  =  -  !i^ 

in  D  and  c^p  must  vanish  on  L.   Along  S  the  conditions  are 
given  by 

(4.20)       -1^  (2<,^^<,^^  -  !^  ^^^  .  ,^  =  C^ 

and 


The  elimination  of  t)   from  (4.20)  and  (4,21)  shows  that  the 
normal  derivative  (j)p  along  S  must  be 
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(4.22)       <t)p„(a,cr,0  - — i-^-^ — 


The  constant  7   can  be  fixed  by  observing  that  If  X   is 
the  outward  normal  derivative  of  %   at  the  boundary  of  D 
then 

^^•^^^    If  ^-^   ^P^O    -"h   (P^p)^^^^P  =  f    P^n  ^^ 
D  L+S 

where  s  denotes  the  arc  length  along  the  boundary  L+S  which 
Is  oriented  counterclockwise  as  shown  In  Fig.  2.   This 
Identity  is  a  condition  for  the  existence  of  a  function  "X 
which  has  a  prescribed  normal  derivative  on  the  boundary  of 
D  and  which  satisfies 


tpp  +^  +  X  ^5  =  K(i:,p) 


or,  what  Is  the  same  equation. 


a(pX_g)      a 


3p   +  §Z   'P^c'  '  PK{C-P) 


In  D.   Such  a  function  does  not  exist  If  the  assigned  normal 
derivative  violates 


j^J"  pK(C,p)dCdp  =  J       pX^  ds 


D  L+S 

The  application  of  (4.23)  to  ^^   gives 
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■Fli<r) 


D 


dCdp 
P 


~^ — JT 


or 


(4.2^)       PJCC) 


D 


P   ^7? 


P7 


o 


7q  -a^gh 


0 


where  Ph  Is  the  breadth  of  the  channel  at  r  =  ah.   Equation 

It 
(4.24)  shows  that  either  F-^{(5r)    =  0  or 


(^.25) 


dCdp 


+ 


D 


P^o 
7o  -  a  gh 


=  0  . 


This  Indicates  a  bifurcation  phenomenon.   The  first  alter- 
native leads  only  to  the  motion  of  an  Irrotatlonal  liquid 
stream  and  hence  It  will  be  Ignored  here.   The  second  alter- 
native gives 


(4.26) 


or 


(4.27) 


7, 


D 


o 


dCdp 


+  P 


D 


iKi^)f-^ 


where 


(4.28) 


D 
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It  should  be  noted  that  if  the  cross  section  of  the 
channel  Is  rectangular  with  depth  h  then  (4.26)  gives 

a^gh  £n   -2L. 
(^•29)  ^o  =   ,   a  ;;' 


and  if 


we  have 


a     ^ 


a-1 


(4.30)  2^a5gh? 


^1 


which  is  the  same  as  (4.1).   This  fulfills  the  expectation 
expressed  at  the  beginning  of  this  section.   Equation  (4.50) 
defines  the  velocity 


^a  _   1  Pgh^ 
^o  -----  p  J-^T- 


which  is  greater  in  absolute  magnitude  than  any  of  the 
velocities  cm/p  associated  with  the  sinusoidal  waves  of  the 
linear  theory  discussed  in  Section  3' 

Equation  (4.26)  defines  the  mean  velocity 
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(4.31) 


u  = 
o 


(f>  dCdp 


D 


D 


D 


dCdp 


P 


dCdp 


a^e-h 


D 


dCdp 
P 


D 


dCdp 
P 


+  1 


When  a  is  sufficiently  large  the  toroidal  channel  can  be 
regarded  as  a  horizontal  channel  of  cross  section  D.   For 
the  limit  of  u  as  a  —>■   oo  we  have 


(4.52) 


lim 
a — >-co 


u  = 
o 


gh  /y  dCdp 

'd 


T 


=  -  /gd 


where  d  is  the  mean  depth  of  the  horizontal  channel.   This 
agrees  with  a  result  in  a  paper  by  the  author  [1],  namely 
that  /gd  is  a  first  order  approximation  to  the  speed  of 
cnoidal  or  solitary  waves  In  a  horizontal  channel  with 
arbitrary  cross  section. 

It  can  now  be  seen  from  equations  (4.19),  (4.22)  and 
the  result  (4.25)  that  the  function  ^^   must  satisfy 


^2pp     p     ^2;; 


2~ 


for  (C,p)  a  point  in  D; 
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%<«'"'?' =  -pji;^^^'^i(°' 


D 

and 


2n 


along   L.   This  Implies  that  we  can  write 


FAa) 
a 


where 

(4.35)     ^2pp  +  ^2£  +  ^2CC  =  ~ 

P  ^ 

for   (C^p)   a  point  in  D; 


(4.3^)         ^2p^^^^^  "  '^^ 


and 


(4-35)         ^2n  "  ° 


along  L.   There  is  no  loss  of  generality  if  we  make  f^     unique 
by  specifying  its  value  at  some  particular  point  say 


Tp^{a,0)    =  0 


From  equation  (4.20)  we  see  that 
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3. 


F^io) 


■""   ■   _.T-i  I  /  ~\  '  T=r7*r*jr 


aF  '  I  a  r 
o 


Fn"oT 

o 


or 


(4.36) 


F-[(a) 


F^(_a)^  ^  ^  ^2 


a|X 


F^Ctjr 


Equation  (4.36)  shows  that   n-,  =  Tln(<^)   is  a  function  of  o 
alone;  but  to  determine  il-,(a)   we  must  proceed  to  the  next  step, 
Tiie  third  order  system  shows  that  *   must  satisfy 


(4.37)     *^    +  -^-^  +  ^^^y 
3pp    P     3CC 


F^'  {o)i^^{p,0        F^(a) 

2 


2  2 
a  p 


P 


F^(a)^2(p'^)^r  ^^^    ^2^^^ 


3  2 
H.a  p 


In  D  ,  and   <J)-,   must  vanish  on   L.   Along  S   the  requirements 
'       3n  o  -1 

are 


(4.58) 


2a 


)'^   j_  p(])   (J) 
la  ^  '^  oa  2a 


4Tln           2  ^3ri^  2ti, 

— ^  4)   (D    +  (D^    — i 1 

a        oa  la    oa  \   2  a 

V  a 


4T12  =  C3 


and 


(4.39)    -- 


a 


la  'la    oa  '2a 


■cF  ll^lo 


*3p  +  ''l*2pp 


The  elimination  of  t]   from  (4.58)  and  (4.39)  gives 
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(4.40)         *^     =  _  F>(^        +  2  +  6[i  +  y^"") 

J>p  o      2pp  ^  ^ 


2.    ^1^1 


ixa" 


+ 


2(l+n)G2Tl'         F^V^g^^l'         M- 


F' 
o 


a 


a      2 


Since   the   normal   derivative      *^        Is   zero   on     L   ,    the 

3n  ' 

existence  condition   (4.23)^  for  71=  '^-^    ,    can  bs  written 


D 


^2 
[^(pct^^)  +  ^_(p*^p)]dCdp  =  a^   <l>^p(a,a,C)dC 


As  we  can  see  from  (4.37)  and  (4.4o)  this  condition  requires 


F' 

o 


.3 


^^,(p,OdCdp  .  T)'^' 


D 


D 


=  -  F 


o 


L^C- 


^2ppd^  +  2P  +  6P;a  +  3PkL^ 


2 


2p(l+M.)aC„iT    F'  ^^2 

+  FT~^-^  -  -3-  /   ^2^r  •  V"i  -  ^^^^2 


Since 


.iw'      P 

D 


the  terms  involving  the  unknown  ¥        cancel  and  we  can  write 
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(4.41) 


^2  rr    V'o^P^OdCdp 


H    /        ^    (a,C)dC    -  — -— 


=   2Pa(l+[x)^C2Ti|(a)    -  a 


D 


nj'  (a) 


>2pp(ct,OdC+2P+6P|x+pM.' 


Tl-]^(a)rij^(a) 


which  is  an  ordinary  nonlinear  differential  equation  for  the 
determination  of  i"!-,  (a)  ,  the  first  approximation  to  the  surface 
wave.   The  coefficient  of   r)'''  (a)   can  be  transformed  into 

r^2            rp    ^^(p,C)dCdp 
m  =  ^  /    Tp    {a,^)di;    -    1' 

^1  D 


%  f       P^2^2n^^  ~  4"  fi    ^2^^p"^P^2p^  +  ^Apf^.)]di:dp 


L+S 


"^  D 


~  If  pii^lr    +  ^IjdCdp 


D 


2C    ^2p' 


which  shows  that  m   is  positive.   The  coefficient  of 

o 


r\-^ia)r]^{o)      is 


m-,  =  -a 


/    ^„^^(a,OdC+2P+6f3!X+3PM-^ 
^=1 


C. 


=  -a  '1  /    [1 ^ ^„..(a,0]dC+2P+6pa+3P:a^f 


I 


J 


-2^^ 


=  -apP+5PM.+5P:a^  -  [Vp^(a,C2)  "  ^2^^°^''^!^^ 
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If  the  curve   L  intersects   p  =  a   orthogonally  we  have 
and 


If  In  addition  to  the  above  definitions  of  m   and  m-,  ,  we 

o        1  ' 

define  m„     as 

m^  =  2Pa(l+(x)^G2 


the  equation  (4.^1)  can  be  written 

(4.42)  m^Ti^'  (a)  -  m^Ti^(r)Tij^(a)  +  xn^T\l^[a) 

The  difference  between  this  equation  and  the  classical  equation 
for  the  determination  of  two-dimensional  solitary  and  cnoldal 
waves  lies  only  in  the  values  of  the  coefficients.   If  we  intro- 
duce the  independent  variable  ")(,_=  aa  and  let  a      tend  to 
infinity  equation  (4.42)  reduces  to  the  equation,  given  by 
Peters  in  [l],  which  determines  solitary  and  cnoldal  waves  in 
a  horizontal  channel  with  arbitrary  cross  section. 

If  m-|   is  negative,  as  it  is  when  L   intersects   p  =  a 
orthogonally,  equation  (4.42)  can  be  written  in  the  form 

(4.43)  AnJ'  (a)  =  -311^  (a)  Til  (a)  +  -K^^{{o) 
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where 


o 


Im-L 


is  regarded  as  known  and 

3m 


2    |m^| 

Is  regarded  as  not  known  since  it  depends  on  Cp.  Equation 
(4,^3)  can  be  solved  by  using  the  following  well  known  pro- 
cedure.  An  Integration  of  (4.43)  gives 


A^i(ct)  =    -  —^  +   A^n;^  +  -2" 


and  if  we  multiply  this  result  by  Tl'(a)  and  integrate  again 
we  find 


,3  ^  .  .2 


A[Ti_^(a)]  =  -T]^  +  A^T^  +  A^T]^  +  A;^ 


or 


(4.45)     A[Tij(a)]"  -  iq-^-^^)  ici^+ri^)  i^^+r]j) 

Since   A   is  positive  and   T-ldj)   must  be  real  when   tj   =  0 
the  product   q-,  q^q^  must  be  positive.   At  least  one  of  the 
q. 's   is  positive  and  we  can  take   q-,  >   0.   The  slope   ^■{(cr) 
must  vanish  at  some  value  of  t]        other  than   r).,  =  q,   for  if 
it  did  not   t]   v/ould  have  to  become  infinite  which  is  physi- 
cally impossible.   If   p  =  a  gives  the  mean  value  of  the  free 


3^' 


^ 


surface  then  we  cannot  have   q_  <  0  and  q^  "*  0   for  if  this 
were  the  case   t]   would  have  to  lie  between  the  extreme  values 
of   q-,  ,  -qp,  -q^  and  thus  be  positive  which  is  again  not 
physically  possible.   It  follows  then  that   q^  >  0  with  q^  ^  0 
and  we  can  take   q^  II  Qp* 

If  we  Introduce  the  new  variable  O   through 

(4.46)         Ti^  =  q^  cos^  O  -  q^  sin^  O 

equation  (4.45)  becomes 


(4.47)         da  =  \-^:~     . ^- 

^'^1+^5      /T-2—2-: 

<■  1-k   sm  O 
where 


If  we  take   a  =  0  when   t]   =  q-.   then  0  =  0  when   a  =  0  and 
(4.48)         a  =  ^'^  ^  ^^ 


%+%     ^o        r~T~2~ 
/  1-k     sm      t 


This   defines    the  amplitude   function 


^i+q^ 
O  -  am  ;   a    /-4s—-     ,    "^ 


This  function  and  the  functions 
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COS  n  =  cos  am  (  a  J    u-s         ,    ^ 


["- 


sin  O  =  sin  am  a  y— in~^  '    ^ 
-   sn  j^a  ,-j^  ,    kj 


2  2     ''     ri^^^       \ 

1  -  k'^sln  O  =  dn  ;  G  /■  j^^  -^  ,  k  ) 


are  called  Jacobl  elliptic  functions  and  in  terms  of  the  en 
function  we  have 


2  f    /qi+q^    /qi+q2 


(4.49)     Ti^(a)  =  -q^  +  (qi+q2)cn-  [a  j-^^  ,/^  J 
This  represents  a  periodic  profile  whose  period  /e"  T  is 


/E^-  I 


IT 


4A     /         dt 


^1+^5  ^o   /,  ,2   .2 
-^       /  1-k   sm 


Tr/2 


^  2  .™^^   /  dt 


•^        /  1-k   sm 


k   sin   t 

Boussinesq  [l8]  was  the  first  to  indicate  the  existence  of  two- 
dimensional  waves  of  the  type  (4.49)  and  the  name  cnoidal  waves, 
obviously  suggested  by  the  notation   en   for  the  Jacobi  elliptic 
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function  and  the  term  "sinusoidal",  was  Introduced  by  Korteweg 
and  De  Vrles  [19] • 
In  addition  to 


k  = 


^  ^1+^3 


It  Is  convenient  to  Introduce  the  quantities 


and 


q  =  £hq. 


£  ^    1  + 


Q.-. 


In  terms  of  these  quantities  we  have 


A. 


— pm^r —  ^  ^1 


Tm-L 


(qg+q^) 


so  that 
(4.50) 


lillL   [3-i(k^+l)] 
^^    6af3(l+M.)^ 


Equation  {^  A3)    takes  the  form 


(4.51)    ^^(/^ 0)  -  -^-[(i-ik^)  +  ^Vcn^riyiA ,  A 


As  a  function  of  9      the  period  of  r^A/^  Q)      Is  given  by 


(4.52) 


T  =  2  '^  f 
\    q^ 


TT 


dt 


■  o 


/i  ,  2    r  2  , 

/  1-k   sm   t 
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Since  we  require  the  mean  depth  of   n-i  (/^  '^)      to  be  zero  we  must 
have 


[ 


T 


^-  o 


/,  „,  2^  ,  .2,  2   2  f'O      ql 
(l-£k  J  +  i  k  en  I  p"  ,  iT 


.k) 


d0  =  0 


and  this  leads  to  the  relation 


(4.53)    (l-«)  ^"^ 


dO 


+  i 


TT 


1-k^  sin^  O  dO  -  0 


o  /l-k^  sin^  Q 


In  terms  of  the  tabulated  complete  elliptic  integrals  of  the 
first  and  second  kinds,  namely 


7r/2 


r/2 


2   .  2 


o    /,  ,  2   .  2 
V  1-k  sm 


r  n  /  d      1 — 

E^(k)  =  /     /1-k"  sin"  t  dt 
'  o 


the   equations    (4.52)    and    (4.53)    are   respectively 


(4.54) 

and 

(4.55) 


T    /qT  _        .    s 

T  <hA  "  ^r^^ 


(l-i)F^(k)    +    I   E^(k) 


0 


The  quantities   k  and  I      can  be  determined  from  these  equations 
when  q  ,    the  height  of  the  crest  of  the  first  order  wave  above 
the  mean  free  surface,  and  the  period  T  are  given. 
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The  depth  of  the  trough  below  the  mean  free  surface  Is 


q  =  ehqp  =  q(-Ck  -l) 


Tnls  shows  that  beside  the  combination   {q,T}   the  combina- 
tions  (q,q]   and   (q,T}   also  suffice  to  determine  r\^. 
The  pressure  equation  (4.4)  can  be  written  as 


Trip, 0,0    =    -p    - 


2p 


2e  I  p     C 


and  by  using  the  foregoing  results  this  Is  the  same  as 


ir{p,a,C,)    =    -p    - 


02 

00 

+ 

^n 

+ 

F. 

2p' 

1 

oa   la 


+  o(e) 


where  the  notation  G(e)  =  0(e)   means    11m 

e  ->  o 


(e) 


0. 


5 .   Summary  and  Discussion 

In  terms  of  the  original  variables  our  theory  up  to  this 
point  shows  the  following.   The  equation  of  the  free  surface 


is 


(■3.1)  r-=a   +   f   =  a   +  hTi=a  +   ehr]-^  +   o(e) 


=  a  +  q 


(1-£r2)    +    tt2on2(|y5J,    ,) 


+   o(q) 
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where  en  denotes  the  Jacob!  elliptic  function.  The  defini- 
tions of  the  parameters  which  appear  in  (5.I)  and  succeeding 
equations  are  listed  below. 

The  clockwise  tangential  speed  is  given  by 


v-L  =  •/g'^  u  -  -/gh  -~- 


r 


-   '"'^^  .  h  ('  *  +£(t>,  1  +  o(e) 
r      I  oa    ia  ' 


'gn 
r 


.  h[F^(a)  +  eFj(a)]  +  o(e) 


^^-  .  h 
r 


eC^a^    eF'(a)ri  (a) 

F'(a)  +  =^7--^ ^ ^— - 

o      F  '  I  a  j       ajx 
o 


+  o(e) 


and  from   (4.27),    (4.50)   and    (4.51)    this    is 


1   +   qU-^l  [3--«(k"fl)] 
h   6a^^Lp(l+M.) 


(5.2)      v^   =    -   |/iv 


^h 


,'l+M- 


\      hap. 


(l-^k2)+ik2cn2C|     a|:,kl 


■^  +   o(q) 


The  radial   speed   is   given  by 


_   /gh   „   _  v/gh         .      _   i/gh    r^2 


V^    =      /r~   V    = 


2. 


2   -/£-"-   /T-   •      p    -  /e 


=i   [£^(t,2^^+o(e^)] 


-   i^  £^/2   _1__  ^      (p,^)    +   o(e2)] 
a  ' 


rv-^i  I  I 


a   LL 
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that  Is, 


(5.3)   V 


/ 


2 


^     / 
a    [I       Y 


en 


k 


■nf^   1^  kV 


•dn 


( I  j^ .  •^';  -  °(^'/^'  • 


The  axial  speed  is  given  by 


v^  = 


3 


/^ 


w 


2, 


^    [b2^^,    +  o(e^)] 


which   is 


(5.4)      .     ,/f-M^-    i£.^(,.,3/2,      (p,,) 


2  WhA 


•  en     -  I    /^   ,  k   i     sn     ^  I    /^    ,  k  )  •      dn  (  - 


2  V  JhA 


e    M  ,k:  +  o{q'/2) 


2  y  hA 


For   the  pressure  we   have         ,    o 

(F    )'^ 

^    o 

-    p    - 


p   =    5gh7r  =    5gh 


2p 


2-   +   ^1 


r                 o    1    1 
+  e[  Cg —    J 


o(e) 


or 
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r 


2 


2   q|mJ[3-^(k^+l) 
a  -  r  +  (1  -  \)  .   ^- 
r 


6ap(l-Hi)' 


(5-5)   P  =  5g   / 


+ 


a 


(1  -  ^) 


y   +o(q) 


+  ^  q[(l-^k2)+^k2cn2(|J||  ,k)] 


For  convenience  we  restate  the  definitions  of  the 
terms  In  the  above  formulas  which  do  not  have  obvious 
meanings : 

q  denotes  the  height  above  the  mean  free  surface, 
r  =  a,  of  a  crest  of  the  first  order  wave. 

k  denotes  a  number  which  Is  related  to  tlie  period  T 
of  the  wave.   This  period  Is  given  by 


TT 


T  = 


dt 


0   /l-k^sln^t 


The  equation 


TT 


TT 


(1-^) 


do 


+  I 


Tl-k^sln^O  do  =  0 


0  vO.-k^sin^O      0 


relates  the  quantity  I    to  k.   It  Insures  that  the  mean 
free  surface  Is  r  =  a .   Th.ls  equation  with  the  equation 
for  the  period  fixes  I   and  k  when  q  and  T  are  prescribed 
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A  is  defined  by  the  ratio 

A  =  


where 


"^0  =  ^  JJp'^'''h   "^  "^gp^  ^^  ^P 


a   p 


and 


?2 


'J  / 

m^  =  -a  p[3+5[i  +  J)\i   -  -   j      i^^^^ict,^)    d^]  . 

D  stands  for  a  dlmenslonless  radial  cross  section 
of  the  equilibrium  position  of  the  liquid.   It  is  bounded 

by  the  segment  S  defined  by  p  =  a,  C.^   <   1^    <   1^^;    and  the 

1   ^3 
channel  wall  curve  L:  —   p  =  a  +  -p^  ^^T7~^  ^  ^  +  y(C)  • 

•^^(p,C)  denotes  the  solution  of  the  system 

*2pp  +  T^  +  *2?C  '  4  (p.OlnD 

^2p(a'^)  =  ctfi  (a,C)  on  S 


iJ^{a,0)    =   0  . 


LL  is  defined  by  the  ratio 

1  rr    dC  dp 
D 
where  p  =  ^2"^1- 


^  =  P  ./J     P 
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We  do  not  give  the  details  here  but  the  following 
statement  can  be  proved.   In  the  limit  as  the  radius  ah 
of  the  free  surface  is  allowed  to  tend  to  infinity  and 
simultaneously  cause  the  wavelength  hxT   to  become  infinite, 
the  formulas  listed  above  reduce  to  those  found  by  the 
author  in  [1]  for  irrotational  solitary  waves  in  a 
horizontal  channel  with  arbitrary  cross  section.   We  note 
that  if  hxT  is  infinite  then  k  =  1  which  implies  qg  =  q    . 
For  k  =  1  we  have 


cos  O  =  en  U  /-^  ,1  =  secha/-j^ 


and  consequently 


,2  a    /  V^ 


T]  (a)  =  -  qg  +  (q-^+qg)  sec,   2   7   A 

If  we  Impose  the  condition  i-j^(oo)  =  0,  as  we  may  without 
loss  of  generality,  then  q^  =   0  from  which  It  follows 


that  q  =  0,  £   =   1,    and 


■n-,(o)  =  q^  sech  -^ 


,„2  a      /^ 


which  gives  the  profile  of  the  solitary  wave. 

When  the  cross  section  of  the  toroidal  channel  is 
a  rectangle  of  breadth  b  =  ph  and  depth  h  the  quantity  p.  is 

H  =  In  (-^) 
"^       a-1 


^5 


and  ^p  must  satisfy 


^2d  ct^     a-1  <  p  <  a. 


2p 


^2p(-'^^  =   ^-^^) 


^2(ct,0)   =  0  . 
For  this  case  It  is  easy  to  see  that 

2 
(S.6)   ^2(p,C)  =  ^2^P^  =^  |-  ^^^^^^^  -  ^""^^A^^  ' 

and  that  the  motion  of  the  liquid  in  the  channel  is 
two-dimensional.   Associated  with  (5-6)  we  have 

a-1 

SO  that  A  is  given  by 

a  ln5(-«^) 
,  _  g-l  

=  5  +  3  in(^)  +  3  In^(^) 
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For  the  case  of  a  rectangular  cross  section  the 


mean  tangential  speed  over  a  period  is 

T  a+f 


where 


Vn 


a-h 


V,  r  dr  dO 


T  a+f 


o  a-h 


r  dr  d9 


^J^ 


1+U- 

1  -  ^ 
^        2a 


<1  + 


[5+'"-H+3kL^][5-^(k^+l)] 
6  atx(l+tJ.) 


}  +o(q) 


p,  =  In  TT^  =  In 


a-1 


a-h 


In  accordance  with  one  of  Stokes  '  definitions  for  the  speed 
of  a  two-dimensional  undulatory  wave,  |  v-,  |  can  be  regarded 
as  the  velocity  of  a  cnoidal  wave  on  a  cylinder.   In  the 


limit  as  a  ->  CO  , 


V-,  I  reduces  to 


lim 
a— >oo 


Vn 


^  (1  +  f^) 


which  Is  the  well  known  formula  for  the  approximate  speed 
of  a  two-dimenslona]  solitary  wave. 

From  the  two-dimensional  theory  of  Section  J>   the 
wave  length  of  the  longest  characteristic  sinusoidal  wave 
is  27ra  and  the  corresponding  free  surface  wave  profile  is 
given  by 

r  =  a  +  A.,co3  9  . 

This  represents  a  cylindrical  surface  whose  cross  section  has 
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a  single  crest  and  a  single  trough,  but  the  height  of 
the  crest  above  the  mean  free  surface  Is  the  same  as  the 
depth  of  the  trough.   For  a  toroidal  channel  with 
rectangular  cross  section  the  onoidal  wave  of  length  2ira. 
defined  by  (5.I)  Is  also  a  cylindrical  wave  with  a  single 
crest  and  a  single  trough  but  for  this  wave  the  height  of 
the  crest  is  not  necessarily  the  same  as  the  depth  of  the 
trough.   Generally  speaking,  the  part  of  the  cnoldal  wave 
above  the  mean  free  surface  is  narrower  in  cross  section 
than  the  part  below.   Furthermore,  if  a  is  sufficiently 
large  the  depth  of  the  trough  of  the  enoldal  wave  can  be 
made  as  small  as  we  please  while  the  height  of  the  crest 
is  kept  finite.   These  observations  can  be  used  to  justify 
the  interpretation  of  the  cnoldal  wave  of  length  27ra  as  a 
two-dimensional  "solitary"  wave  on  a  cylindrical  liquid 
body. 

If  the  cross  section  of  the  toroidal  channel  is 
not  rectangular,  ^pCp^O  depends  on  both  p  and  C  and  since 
each  of  r,  v-,  ,  Vp,  v^  and  p  depend  on  ^p  the  motion  has 
three-dimensional  aspects.   To  the  order  of  approximation 
given  above,  the  height  of  the  free  surface  does  not 
depend  on  ^,  that  is,  the  height  is  constant  in  a  lateral 
direction  across  the  channel.  It  should  be  remarked,  however, 
that  this  is  no  longer  true  for  higher  order  approximations. 
The  next  term  in  the  development  of  the  wave  surface 
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depends   on  r\^  and    (4. 38)    shows   that 


(l-Hi)(l+3(i)Ti^(a)  c^{l-Hi)il+2[i)y]^{a) 


no  = oTTT, + 


2a  p.  a   M- 


(l^H)^c2 


2   +    (n^)c,   -         'cc^{^^^      fUo) 


2  a   M-         '     v-'H-;-3  ^  ^2' 

1         " 

—  n^  (c)  ^2^P'^^ 


a 

The  last  term  in  this  expression  modifies  the  amplitude 
of  T^p  in  the  C  direction.   Therefore  the  approximation 

2         2 
f=  hsT]^  +h  e  T]   +o(e  ) 

furnishes  an  example  of  a  three-dimensional  cnoldal  wave. 

In  conclusion  we  wish  to  point  out  that  although) 
the  foregoing  analysis  is  confined  to  irrotatlonal  motion 
it  can  be  extended  to  study  motions  which  may  exist  as 
bifurcations  from  the  motion  of  a  rotational  liquid  stream. 
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